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Using quantum Monte Carlo (QMC) simulations combined with Maximum Entropy analytic con- 
tinuation as well as analytical methods, we examine the one- and two-particle dynamical properties 
of the Hubbard model on two coupled chains at small doping. The behavior of the single-particle 
spectral weight A(k, uj) as a function of hopping anisotropy t±/t at intermediate interaction strength 
is dominated by the transition from one-band behavior at large t j_ /t to two-band behavior at small 
t±/t, although interaction effects such as band-narrowing, a shift of spectral weight to higher energies 
in the unoccupied antibonding band and reflected structures due to short-range antiferromagnetic 
correlations are also present. A single-particle gap is resolved in the intermediate t±/t Luther- Emery 
phase using Density Matrix Renormalization Group calculations. The dynamical spin and charge 
susceptibilities show features of the expected bonding-band Luttinger liquid behavior, as well as 
higher-energy features due to local excitations between the chains at large t±/t, and evolve towards 
the behavior of two uncoupled chains as t±/t is reduced. For the one hole, large t±/t case, we make 
a detailed comparison between the QMC data and an approximation based on local rung states. 
At isotropic coupling and somewhat larger doping, we find that the dispersion of the single-particle 
bands is essentially unrenormalized from that of the noninteracting system, and that the spin and 
charge response functions have features also seen in random-phase-approximation calculations. 
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I. INTRODUCTION 

In the last five years, there has been much experimen- 
tal and theoretical interest in systems of pairs of strongky, 
correlated-j3pin-l/2 chains.El Materials such as SrC^C^S 
LaCuC>2.5P, and Sri 4 Cu24 42-<5cl contain Cu-0 planes 
which are periodically distorted leading to magnetically 
isolated two-chain CU2O3 ladders with approximately 
isotropic magnetic coupling. These materials are spin- 
liquid insulators, i.e. they have short-range antiferro- 
magnetic order and a gap in the spin excitation spec- 
trum. Recently, the Vanadium-based materials CaV205, 
MgV2C-5 and a'-NaV205 have been proposed to contaip, 
V-0 structures similar to the Cu-0 ladder structures.. 
While CaV2C>5 and MgA^Os are spin liquid insulatorsp 
NaV2C>5 is a one-dimensional antiferromagnet that un- 
dergoes a spin-Peierls transition.! The organic copper 
compound Cu 2 (C5Hi 2 N2)2Cl4 also appears to be a two- 
chain ladder material with antiferromagnetic spin-liquid 
propertiesB Many of the magnetic properties of these 
undoped, insulating materials such as the presence and 
size of the spin gap, the magnetic susceptibility, and the 
magnetic excitation spectrum, are well-desGsjJpcd by the 
Heisenberg model on an isolated ladder .E2Hl2ra However, 
interladder coupling has been argued to be important for 
the MgV2C>5 material)!! and a frustrating next-nearest- 
neighbor interaction as well as a larger interchain than 
intrachain coupling for C^^sH^^^OU-liS 

There have also been attempts to [-synthesize doped 
ladder materials. Hiroi and TakanctJ have succeeded 



in doping the material LaCuC>2.5 by substituting Sr 
for La. While the doped material shows a finite spin 
gap and metallic behavior at low temperatures, there 
is no transition to a superconducting state, which had 
been suggested theoretically.LJlLj However, band struc- 
ture calculations indicate that this material probably 
has a substantial three-dimensional coupling between 
the cha.ins which would make an isolated ladder picture 
invalid. EIlt3 The Sr^Cu^OAg-^ material has been doped 
with Cat3 and Ca and La.cJ At low doping, the spin 
gap present in the insulating parent compound persists 
and the resistivity becomes more metallic. At high dop- 
ing levels and high pressures (3 GPa), the composition 
Sro.4Cai3.6Cu2404j-|84 undergoes a superconducting tran- 
sition at T c = 12K.Q Anomalously, the transition temper- 
ature drops rapidly if the pressure is further increased. 
One issue that remains unclear is the extent to which 
the spin gap persists in the higliporessure superconduct- 
ing composition. Mayaffre et a/.Eil have found low-lying 
spin excitations in the Cu-0 ladders in the normal state 
of the superconducting composition Sr2Cai2Cu2404i us- 
ing 63 Cu-NMR measurements, perhaps indicating two- 
dimensional behavior in the high-pressure superconduct- 
ing phase. 

Although the a'-NaA^Os material was originally 
thought to consist of charge-disproportionated one- 
dimensional V-0 chains, recent X-ray diffraction mea- 
surements and band structure calculations^ suggest that 
it may also be a doped ladder system. The material 
would then have one electron per V-O-V rung, corre- 
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sponding to a quarter-filled conduction band, and an in- 
terchain hopping that is two to three times the intra- 
chain hopping. A quarter-filled two-chain ladder system 
with large interchain hopping can form an effective one- 
dimensional antiferromagnet via the interrung magnetic 
coupling of electrons in the rung bonding state. 

There is growing evidence that a spin-liquid state like 
that found in the ladder materials is closely connected to 
the physics of the high-T c superconductors. For exam- 
ple, the high-T c compounds show signs of a spin gap 
above the Neel temperature and above the supercon- 
ducting transition temperature at small doping.E3 The 
"stripe phases" of the-,two-dimcnsional CuC>2 planes in 
the high-T^ cuprates,cJ which recently have received con- 
siderable attention, are in some sense, ladder-like struc- 
tures. Also, the copper-oxide-based ladder materials are 
clearly closely related to high-T c materials, since the 
ladder-like arrangements of Cu-atoms are produced by 
line defects in CuC>2 planes. 

The success of the two-chain Heisenberg model calcu- 
lations at explaining the magnetic properties of the in- 
sulating materials has lead to the consideration of the 
two-chain Hubbard and t—J models as a starting point 
for the doped materials. In fact, predictions of possi- 
ble properties based on these models predated the ex- 
istence of doped materials. The most interesting of 
these predictions was the possible existence of a spin- 
gapped superconducting state with a d x 2_ y 2-\ihe order 
paramelerJiii isotropic or near-isotropic materials at low 
doping.E3 t3 While the presence of superconductivity in 
the Sro.4Ca13.6Cu24O41.84 material seems to bear out this 
prediction, the extent to which isolated ladder models are 
directly applicable to doped materials is not yet clear. 
The assumption of isolated ladders is probably not justi- 
fied in the doped LaCu02.5 material, and the absence of a 
spin gap in Sro.4Ca13.6Cu24O41.84 does not fit in with the 
picture obtained from the two-chain models. Within the 
two-chain Hubbard model, the relative strength of the 
superconducting correlations, which can at most follow 
a power law decay in one dimension, is quite sensitive to 
the model parameters, especially the doping level and the 
hopping anisotropy.ES Since these parameters also vary 
from material to material, it is crucial to explore fully 
the dependence of a wide variety of properties of isolated 
ladder models on them. Here we will consider the be- 
havior of the single and two-particle dynamical response 
functions of the two-chain Hubbard ladder as a function 
of doping and hopping anisotropy. 

The Hamiltonian for a ladder consisting of two coupled 
chains of length L can be written 

+U^2n t ^n h xi, (1) 
iX 

where A = 1, 2 numbers the chains and i — 1 . . . L the 
rungs of the ladder, t is the hopping matrix element 



within a chain, and t±_ the hopping matrix element be- 
tween the chains. 

There has been much work published in the last few 
years on the behavior of the one-band Hubbard model 
on ladder systems. One important approach is a weak- 
coupling method in which the ground state properties 
of the system are investigated as a function of t± and 
U starting from the limit of weak U, using renormal- 
ization group techniques (see, for example, Refs. 26|,p7|). 
Particularly relevant for realistic two-chain model caL- 
culations is the phase diagram of Balents and FisherEa, 
which is calculated in the limit U — *■ 0+, but for arbi- 
trary values of t± . The resulting phases can be classified 
by the number of gapless spin and charge modes within a 
bosonization picture. Since symmetric and antisymmet- 
ric (fcj_ = and k± = n) modes are possible between the 
chains, there are 9 possible phases, labeled CmSn in Ref. 
p8| , with m = 0,1,2 the number of gapless charge and 
n = 0, 1, 2 the number of gapless spin modes. For exam- 
ple, the CI SO phase has one gapless charge mode and a 
spin gap. This classification scheme has become widely 
used to describe the phases of one-dimensional coupled- 
chain systems in general since it provides a useful com- 
mon notation. Of the nine phases possible in this classi- 
fication scheme, seven are present in the weak-coupling 
phase diagram of Ref. 

One phase of particular interest occurs for most of 
the region of the phase diagram for which two bands 
cross the Fermi energy, i.e. t±/t < 1 + cos7r(n) for 
the infinitely large noninteracting system. The lightly 
to moderately doped isotropic system falls within this 
region. This phase has one gapless charge mode and 
a finite spin gap (C1S0), and is analogous to a one- 
dimensional Luther-Emery liquid. A Luther-Emery liq- 
uid has power-law decaying charge-density-wave corre- 
lations (CDW) and singlet superconducting correlations 
whose relative-strength is governed by a non-universal 
parameter.E3~Eil With dominant superconducting corre- 
lations, the Luther-Emery liquid is a one-dimensional 
analog of a superconductor, with long-range supercon- 
ducting order possible given three-dimensional crossover 
behavior at low temperatures in the real materials. 

Numerical studies of the two-chain Hubbard model us- 
ing the density-matrix renormalization group (DMRG) 
techniqueEJ'LJ show that for a realistic value of the 
Coulomb interaction, U = 8t, the doped system does 
indeed exhibit a spin gap and metallic behavior, i.e. a 
vanishing charge gap, in approximately the same region 
of the t±-(n) phase diagram as the C1S0 phase. in the 
weak coupling calculation of Balents and Fisher £3 In ad- 
dition, there is a direct correlation between the strength 
of the paic,-pair correlation function and the size of the 
spin gapJla The pair order parameter possesses a mod- 
ified d x 2_ y 2 symmetry and the pair correlation decays 
as a power law. These results are consistent with those 
for a t-J model in the small J/t region where-the strong 
U mapping of the Hubbard model applies. lj However, 
the pairing correlation function does not seem to decay 
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slowly enough for isotropic coupling between the chains 
in order for the superconducting correlations to be the 
dominant correlations. 

Dynamical response functions are important both be- 
cause they can be directly probed by experimental tech- 
niques and because they provide a detailed picture of the 
excitation spectra at different energy scales, which can of- 
ten be understood qualitatively or quantitatively through 
various theoretical pictures or approximations. In partic- 
ular, the single particle spectral weight function A(k, ui) 
can be measured directly in photocmission or inverse pho- 
toemission experiments, and the two-particle response 
can be probed using, for example, inelastic neutron- 
scattering experiments. 

Previously, work on the dynamical excitations of the 
doped two-chain system has concentrated mostly on the 
t-J model. In Ref. [33], the photoemission spectra and 
dynamical spin and charge excitations in the half-filled 
and lightly doped t-J ladder were studied using exact 
diagonalization. The starting point of large interchain J 
was used to understand many features of the spectra and 
the isotropic case was studied extensively. In Ref. |3^, 
the single-particle spectral weight of the lightly doped 
isotropic t-J ladder was studied using exact diagonaliza- 
tion, with an emphasis on magnetic features present in 
the single-particle spectrum, and the effect of an inter- 
chain Coulomb repulsion. Recently, Ricra, Poilblanc and 
Dagotto have considered the the effect of the anisotropy 
on the photoemission spectrum in both tlie Hubbard 
a,n(Lt-J models using exact diagonalizationEj Endres et 
alx2 studied both the single particle and two-particle 
dynamical correlation functions of the half-filled Hub- 
bard model using quantum Monte Carlo (QMC) tech- 
niques and a Maximum Entropy analytic continuation 
(MaxEnt), emphasizing the behavior as a function of the 
hopping anisotropy, and making a comparison to various 
analytic calculations. The purpose of the present work is 
to extend the QMC/MaxEnt calculations to the doped 
case. Some of the work is contained in unpublished form 
in Ref. [37|, and a few of the results were publishedj-in 
preliminary form as part of a conference proceedings. E3 

Here we will calculate both the single-particle spectral 
weight function, A(k, uj), and the dynamical two-particle 
spin and charge susceptibilities, Xs(q, w ) and Xc{<\,u). It 
is important to examine both types of correlation func- 
tions in order to make contact with the weak-coupling 
phase diagram in which, in general, well-defined quasi- 
particles do not exist, but the phases can be understood 
in terms of bosonic spin and charge excitations. The 
single-particle spectral weight makes contact with the 
original band structure and its evolution with interac- 
tion, and the two-particle dynamical correlations con- 
tain information on the nature of the low-lying bosonic 
magnetic and charge excitations of the system. The 
single-particle spectral weight function is defined in the 
Lehmann representation as 



A(k,c) = |^e-^(l + e -^)|a|c k , T |/')| 2 
1,1' 

> 8(u - {E lf - Ek)), (2) 

where Z is the partition function, \l > is the exact many- 
body eigenstate with the energy Ei, (3 is the inverse 
temperature l/fc^T and Ck, CT = Si \ c i,Ao-e* k ' Ri ' x /v / 2i. 
Throughout this paper, we will use the notation k = 
(k±,k) for momenta, where the transverse component, 
k±, can take on the values and n, and k denotes the 
component parallel to the chains. We use the grand- 
canonical Quantum Monte Carlo (QMC) algorithm to 
determine the expectation values of the correlation func- 
tions in imaginary time, and invert the spectral theorem 

. r°° e -™ 

G(k,T) = (c kjT (r)4 )T (0))= / A(k,u)du, 

(3) 

using the Maximum Entropy analytic continuation 
method (MaxEnt). Ej In order to achieve high resolu- 
tion, we use a likelihood function which takes the error- 
covariance matrix of the QMC data-and its statistical 
inaccuracy consistently into accounted, an/i. incorporate 
various moments of the spectral weight .eil While the 
fermion sign problem is present, we have chosen system 
parameters, i.e. U /t, filling (n), and lattice size L so that 
low temperatures can be reached with sufficient accuracy. 
The results presented here are based on QMC data with 
good statistics, i.e. averages over 10 5 updates of all the 
Hubbard-Stratonovich variables result in G(k, t)'s with 
absolute errors less than or of the order of 5 x 10~ 4 . 

The two-particle correlation functions are defined sim- 
ilarly in the Lehmann representation, 

Xs , c (q, w) = i£ e-^< (1 - e-^)\ (l\O s ,M | 2 
1,1' 

■ S(uj - {E v - Et)), (4) 

with 0,(q) = E P (4,+q,T c p,T " c p+ q ,i c P,i) and °c(q) = 
S P a c jp+q <t c p o" ^he calculation in imaginary time and 
the analytic continuation are analogous to the calcula- 
tions for A(k, uj). 

Although the QMC results are relatively high resolu- 
tion, some low-energy properties such as single-particle 
and spin gaps are at an energy scale of the order of J — 
At 2 /U, and are not resolvable with QMC/MaxEnt. In 
order to complement the QMC calculations, we will alsa 
present some Density Matrix Renormalization Groupo 
(DMRG) results for low-energy gaps and static correla- 
tion functionsJla The DMRG can provide very high reso- 
lution results for such quantities, but cannot yet provide 
information on the full dynamical properties. 

This paper is organized as follows. In Sec. O, we 
will examine the general evolution of both A(k, uj) and 
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the case of one 
will be treated at 



Xs,c(<l, as a function of t±/t at U/t = 4 and moderate 
doping, (n) — 0.94, i.e. two holes on a 2 x 16 lattice. 
We emphasize making contact with the band transition 
as a function of t±/t as well as the weak-coupling phase 
diagram. DMRG calculations of the single-particle gap 
and the static spin structure factor will be used to re 
solve low-energy features. In Sec. [] 
hole «n) = 0.9875 on a 2 x 16 lattice, 
U/t = 8 and large t±/t = 2.0 Here a direct comparison 
will be made with a local rung approximation (LRA), 
a strong-cou pling , e xpansion based on the states of an 
isolated rungE£Hl3H. We find very good agreement for 
the two-particle as well as the single-particle dynamical 
correlation functions. Finally, we will treat the isotropic 
system, t±/t = 1.0 at a higher doping (n) — 0.875 and 
U/t = 4.0 in Sec. This case is interesting since the 
system is somewhat further from half-filling, and the an- 
tiferromagnetic correlations are weaker. We will com- 
pare the two-particle dynamical correlations with Ran- 
dom Phase Approximation (RPA) calculations, and the 
single-particle spectral weight function with the nonin- 
teracting bands and with slave boson mean-field calcu- 
lations. We find that the dispersion of the bands is es- 
sentially unrenormalized from that of the noninteracting 
system and that the RPA calculations lead to a qual- 
itative understanding of some features of Xs (l) w ) arL d 
good quantitative agreement with most of the features 
X c (q, w). The points in the zero-temperature (n)—t±/t 
phase diagram at which we make calculations are shown 
in Fig. |l| as points in the {7 = phase diagram. 



II. EVOLUTION WITH t±/t 

In this section, we discuss the evolution of the single- 
particle spectral weight A(k, u>), and the dynamic spin 
and charge susceptibilities Xs(q, and % c (q, w), with 
decreasing t±/t at low doping, (n) = 0.94. We choose 
intermediate coupling, U/t = 4, and show results on a 
2 x 16 lattice for which (n) = 0.94 corresponds to a dop- 
ing of two holes from half-filling. The inverse tempera- 
ture, 1/T = P = 8t is low enough so that the system 
is effectively in the ground state for this system size, to 
within the resolution of the QMC calculations and the 
MaxEnt analytic continuation. 

As we will see below, the behavior of A(k, lo) is domi- 
nated by the band transition as a function of t±/t. The 
dispersion of the bands is only moderately renormalized 
by band-narrowing interaction effects, and the one-band 
to two-band transition takes place at tj_/t = 1.6, a signif- 
icantly smaller value than the noninteracting (L — ► oo) 
value of t±/t = 1.98. There are also additional features 
present due to the interactions. 

We display A(k,u) in Figs. |-| for t±/t = 2.2, 1.6, 1.0 
and 0.2. The band transition can be clearly seen. For 
t± = 2.2, only the k± = bonding band lies below 
the Fermi level, and the dispersion of the clearly defined 



peak follows that of the U = band (shown as a dotted 
line). The Fermi wave vector within this band is given 
by kp = w(n) = 157r/16. In the antibonding (k± = tt) 
channel, the spectral weight is in the unoccupied lo > 
portion of the spectrum and lies mostly above the non- 
interacting band. The total spectral weight is small, and 
while the peaks have approximately the dispersion of the 
noninteracting band, they are shifted to higher energy 
and are quite broad. In addition, there are some re- 
flected structures in this channel near k = probably 
due to short-range antiferromagnetic order. 

At t±/t = 1.6 (Fig. ||) the interacting system is just 
below the band transition. The bonding band is al- 
most completely occupied and the antibonding band is 
almost completely unoccupied, leading to k F as and 
k b F ss 157r/16. There is band flattening in both bands 
near the Fermi points, leading to a higher density of 
states in these regions. This high density of states near 
the Fermi level separated by a wave vector Q = (71", tt) was 
argued in Ref. ^5] to strongly enhance ci x 2_2,2— like pairing 
correlations in this regime. According to the weak cou- 
pling phase diagram of Ref. ^8] and the numerical DMRG 
calculations of Ref. [l5|, this phase is a C1S0 phase, a 
Luther-Emery liquid with a spin gap. Structures resem- 
bling the reflected bands present at half-filling are also 
seen, probably reflecting the short-range antiferromag- 
netic order present close to half-filling. 

As t±/t is further reduced, as seen in Fig. |J, the occu- 
pation of the antibonding band increases and that of the 
bonding band decreases. Below the Fermi level (u> < 0), 
the dispersion of the bands follows the noninteracting dis- 
persion well. There is some flattening of the bands at the 
Fermi points and reflected structures can also be seen. 

At tx/t = 0.2, Fig. H the bonding and antibonding 
bands are almost equally occupied, and the distribu- 
tion of spectral weight is very similar in the two bands. 
There is a lightly weighted broad reflected structure in 
the lo > portion of the spectral weight, and a broad 
incoherent background at low energies for lo < in both 
k± channels. At these parameters, 2k F ss n, a condition 
which leads to a C1S2 phase (i.e. gapless spin and charge 
excitations) in the weak-coupling phase diagram of Ref. 



It is interesting to determine to what extent the char- 
acteristics of the low-energy behavior are reflected in 
A(k, lo). Since some of the phases predicted as a function 
of t± are well-understood one-dimensional phases such 
as a Luttinger liquid (C1S1) and Luther-Emery liquid 
(CI SO) phases, characteristics of these phases should be 
present. In principle, a Luttinger liquid phase should lead 
to two separate anomalous peaks in A(k, lo) for k near kp^. 
corresponding to the separated spin and charge modes.O 
The Luther-Emery phase should lead to a single-particle 
gap of the order of the spin gap, as well as an anomalously 
diverging peak corresponding to the charge mode.Eil How- 
ever, resolving the spin-charge separated peaks with nu- 
merical methods has proven to be rather difficult, even 
for the numerically easier case of the one-dimensional 
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Hubbard modela There are no separately resolvable 
spin and charge peaks or resolvable single-particle gaps in 
Figs. §-§. The DMRG, however, can be used to resolve 
gaps at quite small energy scales. 

In Fig. ||, we display the finite-size scaling of the single- 
particle gap 

Ai = m + - A* - = MN + 1) + E Q {N - 1) - 2E (N) 

(5) 

calculated on 2 x 8 to 2 x 32 lattices with open boundary 
conditions. Here [i + — fi" is the jump in the chemical 
potential at the Fermi surface, and E (N) is the ground 
state energy with minimum \Sz\ and N electrons at a 
particular system size, and the overall density n = N/L 
is kept fixed for different system sizes. The lines are 
least-squares fits to quadratic polynomials in 1/L. We 
have chosen somewhat larger doping, (n) = 0.875 than 
the QMC calculations in order to have more system sizes 
available for the finite-size scaling, and U/t = 8. The 
values t±/t = 1.3, 1.6 have been chosen to be just below 
and just above the two-band to one-band transition, as 
determined from DMRG calculations of the spin gap and 
the pair correlations S3 As can be clearly seen, there is 
a finite value of Ai at L — > oo for t±/t = 1.3, in agree- 
ment with the prediction for a Luther-Emery system in 
Ref. [nj. For t±/t = 1.6, Ai scales to zero as L ~> oo, as 
is predicted for the Luttinger liquid (C1S1) phase which 
should be present at these parameters. 

The behavior of the spin and charge modes can 
more directly examined through the dynamical spin 
and charge response Xs(q 5 ui) and Xc(q, We exhibit 
QMC/MaxEnt calculations of these quantities in Figs. 
fj\ for the same parameters as in Figs. 

For tx/t = 2.2, Fig. |^, the system is predicted to be a 
Luttinger liquid (C1S1 phase) in the bonding band since 
the antibonding band is completely unoccupied. There- 
fore, one would expect to see gapless spin and charge 
excitations in the q± = channel at q — 2k b F — 2ir(n). 
The presence of such low- lying excitations at q = n/8 
(where the wave vector has been wrapped back to the 
first Brilloum zone) can be seen in Fig. [j] (a) and (c). 
The dispersion is consistent with a linear increase, as 
would be expected in a Luttinger liquid, although the 
the resolution is not sufficient to resolve different spin 
and charge velocities. The presence of a spin correlation 
at q = (0, 2k F = ir/8) can also be seen as a small peak 
in the q± = branch of the spin structure factor, S(q), 
plotted in Fig. [HJ Here we have calculated 5(q) at T = 
on a 2 x 32 lattice using the DMRG. Note that the much 
more heavily weighted q± — n branch of S'(q) is flat, indi- 
cating very strong antiferromagnetic correlations across 
a rung of the ladder, but very short-range correlations 
along the legs. 

The majority of the spectral weight, however, is asso- 
ciated with high-energy features in the q± = ir channel 
for both Xs(q,w) and Xc(q,u>). For Xs(<l,u), Fig. |(b), 
there is a dispersive, cosine-like band with larger spectral 



weight near q — and a minimum near q = it. This dis- 
persive band is also seen in the half-filled systemL3 and 
can be understood in terms of a local rung approximation 
in which a triplet excited state on a rung is delocalized in 
a cosine band of width J = 4t 2 /U. This picture, valid for 
large t±/t at half-filling, provides a good description of 
the q± = 7r spin excitations because of the proximity of 
the system to half-filling and due to lack of low- lying spin 
excitations in the antibonding channel. The size of the 
q±_ = 7r spin gap is determined by the band separation 
ij_ — t± c where t± c /t 1.6 for these parameters. 

The q± = ir charge susceptibility, Fig. 0(d), is also 
very similar to that of the half-filled system. There is a 
strongly weighted peak at q = at u/t « 7 which can 
be understood in terms of a local charge excitation on a 
rung. As q is increased, this charge excitation becomes 
a broad band with decreasing spectral weight with in- 
creasing q. The peak in intensity has a dispersion which 
increases with increasing q. There is, however, a small 
amount of spectral weight spread over a region with a 
lower bound which decreases with increasing q. This 
broadening of the region of charge excitations can be un- 
derstood principally in terms of particle-hole excitations 
between the single-particle bands, although vertex effects 
do also play a role. A more detailed description of these 
effects for the half-filled case is given in Ref. [$(| 

At t±/t = 1.6, the system is just inside the doped spin- 
liquid (C1S0) phase, and a spin gap is present. As can 
be seen in Fig. || (a), the q± — branch of % s (q, w) has 
gained a little in spectral weight relative to t±/t = 2.2 
and the band of charge excitations has become narrower 
in u. While there is some spectral weight near u = 
and q = tt/8, as was seen for t±/t = 2.2, here one would 
expect a gap in the spin excitations. The finite reso- 
lution here, however, precludes resolving a gap of the 
order of O.lt or less. In the qj_ = n branch, Fig. || (b), 
the cosine-like dispersive band is still present and has 
the same bandwidth, 2t, as for t±/t — 2.2, but is shifted 
downwards in uj so that the q = n peak is close to lu = 0. 
This shift of the band is expected since the gap in the 
q = 7T excitations in the large t± regime is determined by 
the band separation, and should go to zero at the band 
transition. The minimal spin gap, of order O.lt, should 
occur at q = (tt, k F + k F = 157r/16), but the limited 
resolution in q and in to does not allow us to unambigu- 
ously resolve the position of of the minimum, or the size 
of the gap. The spectral weight at q = tt is larger than 
for t x /t = 2.2, Fig. (b). 

The behavior of the spin structure factor S(q), shown 
in Fig. [n], reflects the qualitative change in the spin cor- 
relations due to the band transition. For t±/t = 1.6, 
the peak at q = (0, 7r/8) present for t±/t = 2.2 which 
was due to the 2k F correlations in the C1S1 phase has 
disappeared, as one would expect. There is a larger am- 
plitude in the q± — correlations, reflecting stronger in- 
trachain antiferromagnetic correlations, and a relatively 
broad peak at q = k b F + k F = 157r/16 has developed. 
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This peak reflects short-range incommensurate spin fluc- 
tuations present in the C1S0 phase. El 

The dynamic charge susceptibility, Fig. || (c) and (d), 
has only minimal qualitative changes from the tx/t = 2.2 
case. In the q± = channel, x c (q, to) has a minimum in 
the dispersive band near q = tt/8, with a dispersive band 
which is somewhat broader and more heavily weighted 
than for t±/t = 2.2. In the q± = it branch, the struc- 
tures are similar to the t±/t = 2.2 case, but with some- 
what broader peaks, and, as expected, are also shifted to 
lower lo due to the smaller separation of the bonding and 
antibonding bands. Note that selection rules forbid there 
being any spectral weight at q = (0, 0) for both x c (q, to) 
andxs(q,^)- 

As t±/t is further lowered to 1.0, Fig. [9] and 0.2, Fig. 
|To| , x s (q, lo) and XcClj to) evolve towards a form expected 
for two uncoupled chains. One can see that as t± is low- 
ered, the q±_ = and q± = tt branches of both the spin 
and charge susceptibilities become similar to one another, 
until at t±/t = 0.2, the branches are almost identical to 
within our resolution. At t±/t = 1.0, one sees that the 
dispersion of the peaks in X s (q, w ) nas minima near lo = 
at q = (0,0) and near q = (tt,tt), with heavy spectral 
weight at q = (71", tt). The isotropic system should be in a 
spin-gapped C1S0 phase. The existence of a spin gap has 
been confirmed numerically with DMRG calculationsEj 
for U/t = 8 and is of the order of 0.08i at (n) = 0.9375. 
Again, this spin gap is smaller than the resolution avail- 
able in the QMC/MaxEnt calculations. 

The charge susceptibility, Fig. || (c) and (d), has a 
broad, heavily weighted dispersive band with a minimum 
at q = and lo ~ and a maximum near q — tt in the 
q± = branch. In the q± = tt branch, there is also a dis- 
persive band with a broad peak whose position increases 
monotonically with increasing q. However, the minimum 
position at q = is at lo ~ 3t. In addition, there are 
more lightly weighted dispersive structures at smaller to 
that reach to = at a number of finite wave vectors in 
both q± branches. The overall spectral weight distribu- 
tion and the position of these minima can be understood 
by comparing with Random Phase Approximation (RPA) 
calculations of % c (q, w). Such a comparison will be car- 
ried out in detail for a larger doping, (n) = 0.875 in Sec. 
IV , where we have better resolution data on a larger lat- 
tice. I— I 

At t±/t = 0.2, weak-coupling theorycS predicts that 
the system is close to the line at which the bonding band 
is half-filled at U = 0, and Umklapp processes in the 
bonding band lead to a C1S2 phase, i.e. a vanishing 
of the spin gap. DMRG calculations have shown the 
spin gap vanishes along this line even at large interaction, 
U/t — 8. Figs. 0(a) and (b) are consistent with two van- 
ishing spin modes, to within the resolution of the QMC. 
Both the q± = and the q± = tt branches have a sine-like 
dispersion with minima at q = and near q = it. The 
spectral weight distribution is similar in the two bands, 
except for somewhat larger weight near q = (tt, tt), due 
to the remaining antiferromagnetic correlations, and zero 



weight at q = due to selection rules. The dispersion 
and weight distribution in each— channel are similar to 
that of the ID Hubbard model.t3 The charge and spin 
velocities of the ID Hubbard model, obtained from Ref. 
[46| , are indicated in Fig. [TTJ as lines at low q. As can be 
seen, the \ow-q dispersion agrees well with this ID value 
for both Xs(q 5 to) and x c (q, lo) in both q± channels. 



III. ONE HOLE, LARGE INTERCHAIN 
COUPLING 



For small doping and t± ^S> t, it is possible to ap- 
proximate the properties of the two coupled Hubbard 
chains by starting from a local rung-singlet picture. For 
vanishing doping, the single-particle spectral weight can 
be calculated by considering transitions between a state 
with one hole and a state with two holes. A one-hole 
state with parallel momentum k is given in the LRAl£I 
as \ipi(k)) = Ylt=i e zki \£) /\/L. The expectation value of 
the energy for this state in first order perturbation theory 
is 

E x {k) = (i) X (k)\H\il)i{k)) = (L — l)E a — t± + 2tA cos k, 

(6) 

with A = (l + E b /2t 1 _) 2 /(2 + E^/2t 2 L ). The lowest energy 
state thus occurs for k = tt with the wave vector 



l^ hole ) 



1 



L 

E 



(-1)*\S 1 )\S 2 )...\U)...\S L ), (7) 



and energy El holc = (L - l)E a - t j_ - 2tA. The wave 
function for this state is a product of rung eigenstates, 
all of which have a momentum k± — between the 
chains. Therefore, the total perpendicular momentum 
K± is zero. As discussed in Ref. ^ (see Fig. 7), it is 
energetically most favorable to remove a particle from a 
two-particle rung state \Si) in the k± — channel. The 
removal of a second particle from the one-hole ground 
state will predominantly take place on a rung with two 
particles and have kj_ — 0. (The removal of the particle 
from I le) will also be a fcj^ = process, since the state 
with no particle has k± = 0.) Therefore, we will consider 
only processes in the k± = channel here. The LRA 
state for two holes can then be written as 



|^° le (fc)) 



B 



N 



lholc\ 



where Bn is a normalization constant. Without loss of 
generality, we can choose the direction of the spin so that 
the z component of the total spin vanishes. After some 
calculation, one finds that the dispersion in the k±_ = 
channel for lo < is given by 



.lholc 



(k) = -2tA lholc (l + cosk), 



(9) 



G 



with A lholc = (1 + E b /2t±) 2 /(2 + El/2t\). The factor 
A lholc , which renormalizes the band, is identical to that 
obtained in the half-filled case (i.e. a transition from a 
half- filled LRA state to a one-hole LRA state). In other 
words, to first order in t, the effective bandwidth does 
not change with doping near half-filling. 

In Fig. [l2|, we compare the QMC and LRA results on 
a 2 x 16 system with U/t = 8 and tj_/t = 2.0 for k± = 0, 
part (a), and k± = tt, part (b). The average particle den- 
sity, (n) — 0.96875, corresponds to a doping of one hole 
on a 2 x 16 lattice. The inverse temperature, (3t = 8, is 
the maximum that could be reached due to the fermion 
sign problem, but is nevertheless low enough so that the 
2 x 16 system is effectively in the ground state, to within 
the resolution of the MaxEnt data. 

The solid line in the k± = channel shows the results 
from the LRA calculation. As in the half-filled case, the 
dispersion is in very good agreement with the dispersion 
of the main peak in the QMC data, even at t±/t = 2.0. 
Note that the [7 = band, shown as a dotted line, has a 
larger width and does not fit the peaks as well. There is 
a broad band between uj ~ At and uj w 8t in the k± = tt- 
channel whose bandwidth is the same as that of the band 
in the k± = channel to within error bars. This band 
is displaced upwards significantly from the U = band, 
indicating substantial deviation from the noninteracting 
band structure. Because the resolution of the MaxEnt 
method becomes less accurate for larger uj values, the 
dispersion for k± = it is not as easily discerned as for the 
corresponding band in the k± = channel. While most 
of the spectral weight is contained in these dispersive 
bands, there are indications of a broad incoherent back- 
ground, similar to those found at half-filling, in both k± 
channels. 

In Fig. [L3| we show the charge and spin susceptibili- 
ties for the same parameters as above. Almost all of the 
spectral weight in the spin response function is contained 
in the q± = tt channel. One can understand this result 
by considering the Lehmann representation of X s (q, uj) 
in the LRA picture. Since for q± = the operator O s (q) 
commutes with the LRA Hamiltonian Hq, the matrix ele- 
ments (Z|O s (q)|i') are diagonal in I and V, the susceptibil- 
ity Xs((q± = 0, q), uj) is therefore exactly zero when Hi is 
neglected (i.e. for t± ^> t). In the q± = tt channel, most 
of the spectral weight lies in a band with a cosine-like 
dispersion (Fig. |l^(b)) between uj w 1.9t and uj w 0.81 
A band with identical form.-tp within the error bars is 
found in the half-filled caseES, and can also be identi- 
fied with magnon excitations with the LRA. The spin 
gap in the qj_ = tt channel is the same as the half-filled 
value, AEg alt ^ mling /t = 0.81 ± 0.02 to within the res- 
olution of the DMC/MaxEnt data. The weak-coupling 
phase diagrams predicts a Luttinger-liquid-like phase 
with gapless charge and spin excitations (C1S1) in this 
parameter regime, i.e. no spin gap. There is some indica- 
tion of weakly weighted low-lying spin excitations in the 
q± = channel of the spin susceptibility shown in Fig. 



|T^(a), supporting the weak-coupling phase diagram even 
for the relatively strong coupling U/t = 8. Note that this 
feature is present with a larger weight for (n) = 0.94 and 
U/t — A, as we have seen in Fig. 0(a). 

In Figs. |l3|(c) and (d), we exhibit the charge suscep- 
tibility Xc(Q: uj) for the same parameters. For momen- 
tum q± = tt, shown in Fig. |l](c), there is a broad, 
q-independent excitation at uj « 9t. This excitation can 
be understood in the LRA picture in terms of charge 
excitations on one rung which are also present in the 
half-filled case. Due to the fermion sign problem, which 
becomes worse with doping, it is not possible to obtain 
as high a resolution as for the half-filled system. The 
excitation spectrum for q± = is, however, quite in- 
teresting in that we find a well-defined uj — > 0-Mode for 
q — > with velocity v c /t = 0.59±0.29, even though there 
should be no q± = spectral weight within the LRA 
picture [O c (q± — 0,q) commutes with the unperturbed 
Hamiltonian Ho, analogously to O s (q± = Q,q) for the 
spin excitations] . The presence of low-lying charge exci- 
tations indicates that a metal-insulator transition takes 
place quite close to half-filling. 



IV. ISOTROPIC INTERCHAIN COUPLING AT 
HIGHER DOPING 

In the following, we consider isotropic interchain cou- 
pling t±_ = t, which is relevant for many ladder materi- 
als, at a larger doping, (n) = 0.875. The results are on a 
2 x 64 lattice with U/t — A and at an inverse temperature 
(3t = 5, which is the lowest temperature we could reach 
for these parameters given the size of the average sign. 



A. The single-particle spectral weight 

In Fig. [li], we exhibit the single-particle spectral weight 
A(k, uj) for this system. Two different structures are ap- 
parent in the data. Most of the spectral weight is con- 
tained in a cosine-band-like structure that extends from 
uj w — 2.8t to uj w 1.6i in the k± — channel and from 
uj w — 0.7t to uj rj A At in the k± — tt channel. There is 
also a broad, incoherent background in both k± chan- 
nels for photoemission (uj < 0) as well as for inverse 
photoemission (uj > 0) that contains much less spectral 
weight. A similar incoherent background has also been 
found in stpdies of the one- and two-dimensional Hub- 
bard modelM 

The dispersive structure which contains the majority 
of the spectral weight has the cosine-like form of the non- 
interacting band structure, with a bandwidth that agrees 
well with the U — bandwidth, shown as a solid line, for 
k±_ = 0, Fig. |l4|(a), except for some deviation in the k± 
channel near k = tt. In order to try to understand band 
renormalization effects due to interactions, one can com -, 
pare with slave-boson mean-field theory techniques .e3>c3 
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However, one finds a- narrowing of the noninteracting 
bands of about 10%£3 and a dispersion which does not 
fit the QMC results as well as the dispersion of the non- 
interacting bands. For higher energies, especially near 
k = tt in the k± — tt channel, the agreement between the 
QMC data and the noninteracting dispersion becomes 
worse. One can understand this discrepancy in terms of 
a spin-wave-shake-off effect, which occurs analogously to 
the half-filled case. At half-filling, a self-consistent dia- 
grammatic calculation of the self-energy within a SDW 
state leads to additional contributions to the spectral 
weight, especially in the center of the magnetic Brillouin 
zone, i.e. at k = and k = tt in both k± channels. These 
additional contributions are due to diaaxam classes in- 
volving transverse magnetic fluctuations. EM For the doped 
system, there are still short-range antiferromagnetic spin 
correlations which are reminiscent of spin-density-waves, 
leading to a spin-wave-shake-off effect in the doped sys- 
tem near half-filling. 



B. Charge and Spin Response 

Here we compare the spin and charge excitations of 
the doped isotropic system to calculations within the 
Random Phase Approximation (RPA). This will allow 
us to gain some qualitative understanding of many of 
the features of the charge response especially. The ba- 
sic idea of the RPA is that the bare Coulomb interac- 
tion V® a , (q, u> n ) is replaced by an effective interaction 
Vj^/ A (q, u> n ) in which screening processes present in the 
strongly correlated many-body system are taken into ac- 
count. Analytically summing the ladder diagrams one 
obtains for the spin and charge response function: 



1 ± E7xo(q,u>) 



(10) 



with 

Xo{ 



(q>"») = G o,kK -^)G£, k+q K). (ii) 



k.i/„ 



We show the RPA result for the charge susceptibility, cal- 
culated on a 2 x 64 lattice with U/t = 4, (n) = 0.875, and 
isotropic t± = t in Fig. |l^. Parts (a) and (c) show the 
result of the RPA calculation for T — > for both values 
of q±, where the poles ofy^ PA (q, w), i.e. the zeros of 
the denominator in Eq. ( |l0|) are indicated by diamond 
symbols (O) and the magnitude of the imaginary part of 
xJ? PA (q, w) by the density of the shading. Figs. |l5] (b) 
and (d) show the results of the QMC/MaxEnt calcula- 
tion at fit = 5. The qualitative agreement between the 
RPA calculation and the QMC/MaxEnt results is quite 
good, to within the noise caused by the finite resolution 
and finite temperature of the QMC/MaxEnt method. In 
the q± = branch, there is a well-defined mode that ex- 
tends from u) = at q = to lo rs 5t at q = tt. By taking 



the slope of the maximum weighted portion of this mode 
as q — > 0, we obtain a charge velocity v c /t = 1.74 ± 0.17. 
The band-like structures in the RPA calculation, which 
contain most of the spectral weight, clearly have a sim- 
ilar form to those in the QMC/MaxEnt data. In ad- 
dition, there is also a broad spectrum below the well- 
defined bands, including gapless excitations at wave vec- 
tors qi = 17"7r/32 and q 2 = 257r/32, which are also present 
in the QMC/MaxEnt data, although not so well-defined 
due to the finite resolution of the QMC/MaxEnt and the 
relatively small amount of spectral weight in these modes. 

The q±_ = tt spectra, shown in Figs. [15] (c) and (d), are 
also quite similar for the RPA and QMC/MaxEnt cal- 
culations. In this channel, most of the spectral weight 
for both methods is contained in a well-defined band-like 
structure which extends from oj w 2.5i at q = to u> ~ 6i 
at q = tt. For the RPA calculation, there is a broad spec- 
trum of excitations below this band with gapless excita- 
tions at q 3 = Htt/32 and q A = 28vr/32. The QMC data, 
Fig. |l^(d), again show very similar structures. 

One can understand the source of the gapless excita- 
tions at qi and q 2 in the q± = channel and q 3 and 
<74 in the q± = n channel in terms of the noninteracting 
band structure of the system. In Fig. pi we show the 
bonding (k± = 0) and the antibonding (fcj_ = n) bands 
defined by e(k) for U — with the Fermi wave vectors 
kp 1 and kp 2 within each band indicated. Two gapless 
intraband and two gapless interband charge excitations 
are possible given this band structure. The intraband 
excitations (q± — 0) occur for q 2 = 2tt — 2kp 1 within the 
k± = band and for q\ = 2kp 2 within the k = tt band 
(wrapping back into the reduced Brillouin zone) . The in- 
terband transitions (q± = tt) can occur at 53 = kp l — kp 2 
and qi = kp 1 + kp 2 . At these four wave- vectors, gapless 
particle-particle excitations can occur, leading to uj = 
poles inv^ PA (q,w) and thus the gapless excitations seen 
in Fig. [l5| (a) and (c). 

In Fig. n7], we show the QMC results as well as the 
RPA calculation for the spin susceptibility Xs{<l,w) for 
the same parameter values as the charge susceptibility 
discussed above, with the results plotted in a similar 
manner. While some of the features of the RPA cal- 
culation are present, here the qualitative agreement with 
the QMC results is not as good. There is a single gapless 
(u> — ► 0) mode present in q± = channel in the QMC 
data that increases with increasing q, reaching u> ~ 1.5t 
at q = tt. The RPA calculation, in contrast, shows 
a broad spectrum which extends to u> w At at q = tt 
with the largest spectral weight contained in two disper- 
sive bands which become gapless at 95 w 137r/32 and 
<?6 ~ 277r/32. The deviation of the RPA results from 
the QMC calculation is probably due to the following: 
first, the relatively high temperature of the QMC sim- 
ulation (fit — 5) leads to a thermal broadening of the 
structures so that individual poles could no longer be re- 
solved, even with exact (i.e. with zero statistical errcp). 
data. Secondly, it is known from DMRG calculations^ 
that there is a spin gap in the isotropic system even for 
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moderate doping. Since the RPA calculation results in a 
Fermi liquid state, there is no spin gap in the RPA re- 
sults, whereas in the QMC data, the presence of a spin 
gap, expected from the DMRG results, could lead to the 
vanishing of spectral weight at q§ and qg . 

In the q± = tt channel, shown in Fig. [I?] (c) and 
(d), the RPA calculation shows gapless excitations at 
q 7 « Htt/32, q 8 w 20tt/32 and q 9 w 28tt/32. In the 
QMC data, one can only clearly see the upper part of 
the dispersive band that becomes gapless at qi and has 
a value u> » t for q — ► 0. Most of the spectral weight 
is, in contrast to the RPA calculation, contained in the 
qx = 7r channel and is concentrated at u> ~ 0.5i between 
q = 37r/4 and q = tt. The overall dispersion of this struc- 
ture is somewhat cosine-like, reminiscent of that seen 
in Fig. |l^(b) for small doping and tj_/t = 2 and thus 
could be be partially due to the propagation along the 
ladder of a rung spin-triplet excitations in a spin-singlet 
background, as described previously in the LRA picture. 
Therefore, the QMC data possibly show some hints of a 
spin gap in Xsi^l^) m both q^ channels. 

V. CONCLUSION 

We have discussed the dynamic single-particle, as well 
as two-particle (i.e. charge and spin) response of a doped 
two-chain Hubbard ladder, calculated using QMC and 
the Maximum Entropy method and concentrating on the 
regime of small doping. The general evolution of the 
single-particle as well as the charge and spin response 
has been investigated for a filling of (n) — 0.94, corre- 
sponding to two holes on a 2 x 16 lattice at U/t = 4.0, 
an interaction strength equal to the one-dimensional 
bandwidth. The behavior of the single-particle spec- 
tral weight A(k, u>) is dominated by a one-band to two- 
band transition as t±/t is decreased from the highly 
anisotropic case. There are additional interaction effects 
such as a flattening of the dispersion near the Fermi level, 
band narrowing, and reflected bands due to antiferromag- 
netic correlations which are especially prominent near 
the band transition. The presence of a single-particle 
gap when both the bonding and the antibonding bands 
cross the Fermi surface, while not resolvable within the 
QMC/MaxEnt data, was confirmed using DMRG calcu- 
lations. 

At large t±/t the dynamical spin and charge suscepti- 
bilities show heavily weighted high energy features cor- 
responding to local magnetic and charge excitations on a 
rung, as well as more lightly weighted gapless modes cor- 
responding to Luttinger liquid behavior in the bonding 
channel. The presence of the gapless spin mode is also 
found in the static magnetic structure factor calculated 
with the DMRG. As t±/t is lowered, the spin response 
evolves first to reflect the strong short-range antiferro- 
magnetic correlations at q = (tt, tt), and then evolves to 



the form of uncoupled one-dimensional chains with simi- 
lar dispersion and weight in both transverse channels for 
small t±/t. The charge response is consistent with a gap- 
less mode in the q± = channel as t±/t is reduced. In 
the q± = tt channel, the charge response evolves from a 
form marked by a high-energy local rung excitation to 
a form with a low-energy mode with linear dispersion 
with a weight distribution similar to that of the q± = 
channel as t±/t is reduced. 

For the case of one hole at strong interchain coupling, 
t± ^ 2t and U/t = 8, the results are well-described by a 
Local Rung Approximation (LRA) , in which the states of 
one rung are treated exactly and the interaction between 
the rungs is treated perturbatively. Both the dispersion 
and the weight distribution in the effective band struc- 
ture, given by the single-particle spectral weight A(k, uj), 
is quantitatively well-fit by a LRA calculation in which 
one or two rung singlets are replaced by states with one 
hole, and then delocalized into a Bloch band. The q± = tt 
dynamic charge and spin susceptibilities are also well- 
described by an LRA picture in which local rung charge 
or triplet excitations are delocalized into a band. How- 
ever, there is a well-defined gapless charge mode in the 
q±_ = channel of the charge response signifying metallic 
behavior, while spectral weight for q± = is forbidden 
by selection rules in the LRA picture. 

For the isotropic system at moderate doping, (n) — 
0.875, we still find the overall quasi-particle dispersion es- 
sentially unchanged. However, a broad incoherent back- 
ground with a width of several t, that is analogous to sim- 
ilar features observed in the one- and two- dimensional 
Hubbard models, is found. The dispersion of the quasi- 
particle-like structure agrees well with that of the U = 
band structure. The general features of the charge sus- 
ceptibility, including the origin of wave vectors at which 
gapless excitations occur, can be understood in terms of 
a Random Phase Approximation (RPA) picture. In con- 
trast, the spin susceptibility shows the signature of short- 
range antiferromagnetic ordering not present within the 
RPA picture. 
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FIGURE CAPTIONS 



FIG. 1. The U = phase diagram in the tx-(n)-plane. In the shaded region, both the bonding and antibonding bands 
intersect the Fermi level, while in the unshaded region only the bonding band is occupied. On the dashed line, the bonding 
band is half-filled. The solid points signify parameter values at which we calculate the dynamical properties in this work. The 
open circles represent points at which we calculate the quasiparticle gap with the DMRG. 

FIG. 2. The single-particle spectral weight A(k,u>) for t±/t = 2.2, U/t = 4 and (it — 8 on a 2 x 16 lattice at a filling of 
(n) = 0.94. The two k± branches are overlayed on the same plot since they are well-separated in energy, with the open and 
filled symbols marking the position of the k± — and k± = tt peaks, respectively, and the error bars indicating the width of 
the peaks. The density of shading indicates the amount of the spectral weight. 



FIG. 3. j4(k,w) plotted as in Fig. 2, with the same parameters, except with t±/t = 1.6. 



FIG. 4. A(k,uj) plotted as in Figs. 2-3, except with t±/t = 1.0. 



FIG. 5. A(k,u}) plotted as in Figs. 2-4, except with t±/t = 0.2 and (a) k± = and (b) k± = tt plotted separately. 

FIG. 6. Finite size extrapolation of the quasiparticle gap, calculated using the DMRG for U = 4, (n) = 0.9375, and 
t±/t = 1.3, 1.6. 



FIG. 7. (a) and (b), the spin susceptibility Xs(qjW) and (c) and (d), the charge susceptibility Xc(q,w) for t±/t = 2.2, 
U/t = 4 and /3t = 8 on a 2 x 16 lattice at a filling of (n) = 0.94. 



FIG. 8. The spin and charge susceptibilities plotted as in Fig. 7, except for t±/t — 1.6. 



FIG. 9. The spin and charge susceptibilities plotted as in Figs. 7-8, except for t±/t = 1.0. 

FIG. 10. The spin and charge susceptibilities plotted as in Fig. 7-9, except for t±/t = 0.2. The straight lines indicate the 
spin- and charge velocities of the corresponding ID Hubbard chain (from Ref. 46). 



FIG. 11. The spin structure factor S(q) for t±/t = 1.6, 2.2 calculated with the DMRG with U/t = 4 and (n) = 0.9375. 

FIG. 12. The single-particle spectral weight j4(k,a>) for t±/t = 2, U/t — 8 and (it = 8 on a 2 x 16 lattice at a filling 
of (n) = 0.96875 (1 hole). The error bars denote the width of the peak obtained using the Maximum Entropy technique. 
The dotted lines indicate the dispersion of the noninteracting, U = bands, and the solid line in (a) indicates the dispersion 
obtained from the LRA calculation. 

FIG. 13. The spin and charge susceptibilities (xs(q,w): parts (a) and (b), Xc(q,to): parts (c) and (d)) for t±/t = 2, U/t — 8 
and fit — 8 on a 2 x 16 lattice at a filling of (n) = 0.96875 (1 hole). The q± — channel is shown in (a) and (c), while the 
q± = 7r channel is shown in (b) and (d) . 
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FIG. 14. The single-particle spectral weight A(k,ui) for t±/t = 1, U/t = 4 and /?f = 5 on a 2 x 64 lattice at a filling of 
(n) = 0.8776 ± 0.0001. The k± = and k± — 7v data are shown in parts (a) and (b), respectively. The solid line represents the 
dispersion of the noninteracting (U = 0) bands. 

FIG. 15. The charge susceptibility Xc(q, for isotropic t± = t with U/t = 4 on a 2 x 64 lattice at a filling of (n) = 0.875. 
Parts (a) and (c) show the results of the RPA calculation for both q± channels at T = (grey-shaded data and diamond O 
symbols), while parts (b) and (d) show the QMC results for fit = 5. 



FIG. 16. Band structure of the interaction free case (U = 0) of a 2 x 64 system with t±/t = 1.0 and a particle density of 
(n) = 0.875. 

FIG. 17. The spin susceptibility Xs(q, for isotropic t± = t with U/t = 4 on a 2 x 64 lattice at a filling of (n) — 0.875. 
Parts (a) and (c) show the results of the RPA calculation for both q± channels at T = (grey-shaded data and diamond O 
symbols), while parts (b) and (d) show the QMC results for (it = 5. 
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